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This paper gives a pedagogic review of the envelope formalism for ex- 
citation of zonal flows by nonlinear interactions of plasma drift waves 
or Rossby waves, described equivalently by the Hasegawa-Mima (HM) 
equation or the quasigeostrophic barotropic potential vorticity equation, 
respectively. In the plasma case a modified form of the HM equation, 
which takes into account suppression of the magnetic-surface-averaged 
electron density response by a small amount of rotational transform, is 
also analyzed. Excitation of zonal mean flow by a modulated wave train 
is particularly strong in the modified HM case. A local dispersion relation 
for a coherent wave train is calculated by linearizing about a background 
mean flow and used to find the nonlinear frequency shift by inserting the 
nonlinearly excited mean flow. Using the generic nonlinear Schrodinger 
equation about a uniform carrier wave, the criterion for instability of 
small modulations of the wave train is found, as is the maximum growth 
rate and phase velocity of the modulations and zonal flows, in both the 
modified and unmodified cases. 



1. Introduction 

m 

As of January 1, 2006, Wikipedia introduces the term "zonal flow" thus: 

Fluid flow is often decomposed into mean and deviation from the 
mean, where the averaging can be done in either space or time, thus 
the mean flow is the field of means for all individual grid points. In 
the atmospheric sciences, the mean flow is taken to be the purely 
zonal flow of the atmosphere which is driven by the temperature 
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2 R.L. Dewar & R.F. Abdullatif 

contrast between equator and the poles. 

In geography, geophysics, and meteorology, zonal usually means 
'along a latitude circle', i.e. 'in the east-west direction'. In atmo- 
spheric sciences the zonal coordinate is denoted by x, and the zonal 
wind speed by u. 




Fig. 1. Left panel shows NASA image PIA04866: Cassini Jupiter Portrait, a mosaic of 27 
images taken in December 2000 by the Cassini spacecraft. Right panel shows a simulation 
of plasma potential fluctuations in a tokamak cut at a fixed toroidal angle as produced 
by the GYRO code (courtesy Jeff Candy http://fusion.gat.com/theory/pmp/ |. Note 
that, in the plasma case, zonal flows are in the ^/-direction when slab geometry is used. 



The above definition omits the further qualification that zonal flows 
are restricted to bands (zones) of latitude. This is most clearly seen in the 
banded cloud patterns on Jupiter (see left panel of Fig.QJ where the magni- 
tude, and even sign, of the zonal flows varies with latitude in a quasiperiodic 
fashion. (On Earth, topographic variations like mountain ranges disrupt the 
zonal symmetry.) 

The term "zonal flow" has also recently come to be much used in toroidal 
magnetic confinement plasma physics (see e.g. the review of Diamond et 
alt**) to refer to a mean poloidal flow with strong variation in minor radius. 
The sheared nature of this flow is thought to have the strongly beneficial ef- 
fect of reducing radial transport by suppressing turbulence, thus improving 
the confinement of heat required to achieve fusion conditions. 

The use of the same phrase "zonal flow" in the context of both geo- 
physics and magnetic plasma confinement is no coincidence, as the exis- 
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Zonal flow generation by modulational instability 3 

tence of strong analogies between these fields has become well recognized 
In this paper we work in the plasma context, but point out the relation to 
the geophysical context when appropriate. 

The right panel of Fig.^depicts a section of a simulated tokamak, show- 
ing turbulence excited by gradients in temperature and density (the plasma 
being hotter and denser in the middle section, which has been cut out to aid 
in simulation and visualization). The radial coherence length of the drift- 
wave eddies has been reduced by spontaneously excited zonal flows. The 
magnetic field is predominantly in the toroidal direction, but there is some 
component in the poloidal direction so that the magnetic field lines wind 
helically around the torus, mapping out nested toroidal magnetic surfaces 
that confine the magnetic field in topologically toroidal magnetic flux tubes. 
The helical nature of the magnetic field lines can be seen in the figure from 
the fact that the turbulent eddies have their cores essentially aligned with 
the magnetic field, making the turbulence quasi-two-dimensional despite 
the manifestly three-dimensional nature of the tokamak. 

A generalized polar representation would clearly be most appropriate 
for representing the cross-sectional plane of the torus, but, for the purpose 
of gaining physical insight with a minimum of formalism in this paper we 
use slab geometry. That is, the toroidal magnetic surfaces are imagined 
as flattened into planes, so that Cartesian coordinates, x, y, z, can be 
used, with y and z replacing the poloidal and toroidal angles, respectively, 
and x the minor radius. The x and y directions are indicated in the right 
panel of Fig. ^ The slab approximation is the analogue of the /3-plane 
approximation in geophysics, but note the axis convention is opposite to 
that used in geophysics, with y now the zonal direction. 

Modulational instability of drift waves (the analogue of planetary 
Rossby waves) is a strong candidal El El El El lifer generating these zonal 
flows through a feedback mechanism, in which modulations of the wave en- 
velope excite zonal flows through a nonlinear mechanism (Reynolds stress) 
and the zonal flows enhance the modulation through a self-focusing mech- 
anism. 

It is the aim of this paper to elucidate this theory in a pedagogic 
way using as simple a plasma descrip tion as possible, namely the one- 
field Hasegawa-Mima equation.EII El This provides a simple theoreti- 
cal starting poi nt fo r describing the nonlinear interaction of drift waves 
and zonal flowsPlEl The same equation also describes Rossby wave 
turbulence in planetary flows in the quasigeostrophic and barotropic 
approximations^! HI To emphasize its geophysical connections, we 
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4 R.L. Dewar & R.F. Abdullatif 

shall follow a common practice in the plasma physics literature and call 
the original form of the Hasegawa-Mima equation the Charney-Hasegawa- 
Mima (CHM) equation (although in the geophysical literature the equation 
is called the "quasigeostrophic barotropic potential vorticity equation'^} . 

Unfortunately for the interdisciplinarity afforded by the use of the CHM 
equation, it became recognized in the early '90sPI ESI that a corrected form 
for toroidal plasma applications should be used, which we shall call the Mod- 
ified Hasegawa-Mima (MHM) equation. Although the modification seems 
at first glance to be minor, we shall show that it makes a profound difference 
to the modulational stability analysis because it enhances the generation 
of zonal flows. 

Some early works on modulational instability of drift waves can be found 
in Rcfs. 1171 1181 1191 but these predate the recognition of the need to use the 
MHM equation to enhance the nonlinear effect of zonal flows in a toroidal 
plasma. Both MajumdaJi^ and Shivamoggff^ add a scalar nonlinearity, 
arising from polarization drift and/or temperature gradient, to the CHM 
equation in order to find a nonlinear frequency shift. In the Mima and 
paper, the nonlinear frequency shift comes from time-averaged flow 
and density profile flattening. 

In Sec. (0) we introduce the CHM and MHM equations and in Sec. © 
we introduce the generic form of the nonlinear Schrodinger equation, which 
describes the time evolution of modulations on a carrier wave, and use it 
to derive a criterion for modulational instability. In Sec. JSJ we use the 
MHM and CHM equations to derive the nonlinear frequency shift of a 
finite-amplitude drift/Rossby wave and use it to determine the criteria for 
modulational instability of drift and Rossby waves, respectively. Section © 
contains conclusions and directions for further work. 



2. The CHM and MHM equations 

The Charney-Hasegawa- Mima equation (CHmJHIEI is an equation for 
the evolution in time, t, of the electrostatic potential (p(x, y, i) (or, in the 
Ross by w ave application, the deviation of the atmospheric depth from the 
mearl^J). Here x and y are Cartesian coordinates describing position in a 
two-dimensional domain P, representing a cross section of a toroidal plasma 
with a strong magnetic field, B, predominantly in the z-direction (unit 
vector z). 

In the slab model we take V to be a rectangle with sides of length L x 
and L y . A circular domain would clearly be more realistic geometrically 
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because it has a unique central point, representing the magnetic axis, but 
it is unlikely to add any qualitatively new physics. In strongly shaped toka- 
maks, like the one depicted in Fig. Q one might be tempted to give T> the 
noncircular shape of the plasma edge to add yet more realism. However, we 
caution against this line of thinking because each point in T> represents an 
extended section of field line, over which the drift- wave amplitude is signifi- 
cant. That is, T> does not represent any given cross section of the tokamak, 
but rather a two-dimensional surface in a field-line coordinate space (see 
e.g. Ref . I2U|) . onto which behaviour in the third dimension is projected. 

We assume the ion temperature to be negligible with respect to the 
electron temperature T Cl assumed constant throughout the plasma. The 
strong magnetic field allows the plasma to support a cross-field gradient 
in the time-averaged electron number density, n, and the wave dynamics is 
taken to be sufficiently slow that, along the field lines, the electrons respond 
adiabatically to fluctuations in tp. That is, on a given field line they remain 
in local thermodynamic equilibrium, with distribution function /(r, v, t) = 
constexp(— S/T e ), where T c is the electron temperature in energy units 
(eV) and £ is the total electron energy ^mv 2 — eip, with in the electron 
mass and e the electronic charge. 

Following Hasegawa and Mima, the shear in the magnetic field is as- 
sumed very weak, so that z-derivatives and the parallel component, k\\, 
of the wave vector k can be neglected, thus reducing the problem to a 
two-dimensional one. However, the existence of magnetic shear is crucial in 
one qualitative respect — the foliation of the magnetic field lines into nested 
toroidal magnetic surfaces (x = const in slab geometry). Field lines cover 
almost all magnetic surfaces ergodically, so the constant in the expression 
for the distribution function is a surface quantity. Integrating over velocity 
we find 



where we have decomposed ip into a surface-averaged part, Tp{x, t) = 
P tp(x,y,t) (absorbed into no), and the surface-varying part, tp = Pp = 
<p—Tp. Here we have used the magnetic-surface-averaging operator P defined 
in slab geometry by 



and its complementary projector P = 1 — P. (Note that P and P com- 
mute with d t and V.) Equation Q can also be derived purely from fluid 




(1) 




(2) 
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6 R.L. Dewar & R.F. Abdullatif 

equations, without introducing the distribution function explicitly. 

One can show, by surface-averaging the continuity equation for the elec- 
tron fluid in the absence of sources or sinks, that the surface-averaged elec- 
tron density is independent of t. Thus, to 0((p), hq is independent of t and 
equals the prescribed average density n. This would not be the case if we 
had not subtracted off Tp in Eq. (JIJ , and in this we differ from Haseg awa and 
Mima but follow most modern practice since Dorland et al\^ ED pointed 
out the importance of modifying the electron response in this way. Use of 
Eq. Q leads to what we shall call the Modified Has eg aw a- Mima equation 
(MHM equation). 

Defining a switch parameter s such that s = selects the original CHM 
equation and s = 1 the MHM equation, and a stream function if) = if/Bo, 
we write Eq. (10) of Hasegawa et alU^ as 

£(l*^ + ^-^(* + W)=0, (3) 
dt \ n w ci r e J 

where d/dt = dt + Vg-V, with 

V(fXz . 

v s = =zxVV>, (4) 

i>0 

being the ExB velocity (SI units), ( = Z'Vxv^ = V 2 if> the vorticity, 
u> c i = eBo/rrii the ion cyclotron frequency, and V = +yd y the perpen- 
dicular gradient. As shown in the Appendix of Meiss and HortorP21 this 
is an approximate form of Ertel's theorem for the conservation of poten- 
tial vorticity under Lagrangian advection at the ExB velocity. Note that 
the MHM equation satisfies the expected 21 Galilean invariance under boosts 
in the poloidal direction, y' — y — Vt, E' = E + VB x (so if)' = ifj — Vx), 
whereas the original CHM equation does not and is therefore unsatisfactory 
for plasma physics purposes. 

We now rewrite Eq. @ in a more explicit way [cf. Eq. (1) of Smolyakov 
eta^ 

(d t + v E -V + v„-V)($ + S s ^) - (d t + ve-V)p 2 s V 2 iI> = , (5) 

where the characteristic drift-wave scale length p s = ixT.^X'c/TOi) 1 / 2 is the 
sound speed divided by ui c i, and the electron diamagnetic drift is defined 

a Even in the absence of topography, we do not expect Galilean invariance in geophysical 
application of the CHM equation, as the /3-plane is not an inertial frame. In the plasma 
confinement application, a poloidal boost in polar coordinates would also be to a rotating 
frame, but the slab approximation implies we should ignore any Coriolis effects and 
Galilean invariance in the poloidal direction should apply. 

h v t = v* is the analogue of /3 in the geophysical application of the CHM equation. 
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by 



v* = — 



T c z xV?io 

eB n 



(6) 



The ordering in Ref. 11 makes all terms in Eq. J^J of the same order. Thus, 
balancing dt<p and dtp 2 V 2 ip 2 we see that p s is indeed the characteristic 
scale length for spatial fluctations. Balancing dtf and v*-V<y9 we see that 
the characteristic time scale is p s /v*, and balancing ve and v* we see that 
the characteristic amplitude of potential fluctuations is (T e /e)p s / L n , where 
L n is the scale length for radial variation of no- We assume p s /L n <C 1, 
so the waves have small amplitudes compared with the thermal potential. 
However, with k a typical fluctuation wavelength, and £ a typical dis- 
placement of a fluid element by the waves, can be order unity, and thus the 
equation can describe strong turbulence. 

Projecting Eq. JSJ with P and P we can split it into two equations, one 
for the surface-varying part and one for the zonal-flow part 



In the MHM case, s = 1, Eq. © reduces to Eq. (2) of Ref. d t V 2 ip = 
-Pzx Vf VV 2 f 

Although physically an approximation, we shall in this paper regard 
the CHM/MHM equations as given and treat them as exact equations even 
for the mean flow component of ip, which we assume to vary on longer 
and slower length and time scales than assumed in the maximal balance 
ordering discussed above. 

3. Waves and mean flow 

Assuming there is a scale separation between fluctuations and mean flow, 
we introduce an averaging operation (■) which filters out the fluctuating, 
wavelike component of whatever it acts on, leaving only a slowly varying 
component related to the mean flow. This operation can be realized ex- 
plicitly by convolution with a smooth, bell-shaped kernel of width (in time 
and space) long compared with the fluctuation scale but short compared 
with the mean flow scale. Alternatively we can define it implicitly via the 
test-function formalism introduced in Appendix A of Ref. 22 Either way, 



(d t +zxVfV)(l - p s 2 V 2 )Vi + [v* - zxV((J s , - p s 2 V 2 )#W 
= p 2 PzxV?/>-VV 2 -0 , 



(7) 



d t (S Si0 - pi V 2 )^ = p 2 s Pzx V^-VV 2 ^ . 



(8) 
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Fig. 2. Visualization showing disruption (lower right) of a coherent wave train (lower 
left) through generation of zonal flows by modulational instability in a tokamak simula- 



averaging can be shown to commute with d t and V to all orders in e, the 
ratio of fluctuation scales to mean-flow scales. 

We then split tp into a slowly varying mean flow part, tpo = (V^i an d a 
fluctuating part, ip± = ip — {tp). 

Note that, except when the mean flows are purely zonal, (■) is distinct 
from the surface averaging operation P- we used to set up the MHM equa- 
tion. In this we differ from Champeaux and Diamond*^, who, in effect, take 
P to be the same as (•} irrespective of the direction of the mean flows. 

We consider the case of modulations carried on a coherent wave (see 
e.g. Fig.[2J), rather than a broad turbulent spectrum.EH (As the CHM and 
MHM equations include no drift-wave instability mechanism, the origin 
of this wave is outside the theory — it is an initial condition.) Taking, for 
simplicity, v„ to be a global constant we assume the carrier wave (also called 
the pump wave in some analyses) to be a plane wave and write 



where A is a slowly varying complex amplitude and c.c. denotes complex 
conjugate. 

To begin, we take A and the mean flow, (ve), to be constant and treat 
the carrier wave using linear theory. (Nonlinear effects will be discussed in 
Sec.0]) Linearizing Eq. 10) we find the dispersion relation in the CHM case 




= A(r, t) exp(ik-r — w^t) + c.c. , 



(9) 
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Zonal flow generation by modulational instability 9 

to be 

- = tt^ + tt|^<->' < 10 > 

whereas, in the MHM case, s = 1, and assuming purely zonal mean flow 
((ve) = v E f), it is 

"* = Tri? +k -<^- (11 > 

In the latter case, the mean flow causes a simple Doppler shift of frequency, 
but for the unmodified CHM equation the Doppler shift is reduced by a 
factor p 2 k 2 /(l + p 2 k 2 ). 

We shall use the frequency shift due to mean flow to calculate the 
nonlinear frequency shift. Otherwise we can ignore it. The group veloc- 
ity, v g = divk/dk, in the absence of a mean flow, is the same in both cases 



1 du k y 2p 2 s kk 



y 



u* 5k l + p 2 k 2 (l + p 2 k 2 ) 2 ■ 
We shall also need the dispersion dyadic VkVk^fe 



(12) 



1 d 2 uj k _ Q 4 k y kk 2 ky + yk + fc y l 



Vvdkdk ^ s (l + p 2 s k 2 ) 3 rs (l + p 2 s k 2 ) 2 
where I is the unit dyadic. 

4. Nonlinear Schrodinger equation and modulational 
instability 

We largely follow the sim ple introduction to modulational instability theory 
given in Dewar et al. starting with the nonlinear Schrodinger equation 

| + £.v)A-*,MM-i^.WA, (14, 

where d/dh denotes the gradient in k-space and Alo the nonlinear frequency 
shift, a nonlinear functional of the amplitude \A\ (cf. e.g. Ref. 125(1. If the 
scale length of the modulations is 0(e _1 ) compared with the wavelength 
of the carrier, then the v g -V term on the LHS of Eq. ((141) is 0(e) whereas 
the VkV g :W term on the RHS is smaller, 0(e 2 ). Assuming the nonlinear 
frequency shift to be of the same order, we see that Eq. ((14(1 expresses the 
fact that, on a short time scale, modulations simply advect with the group 
velocity, while on a longer timescale the nonlinear frequency shift causes 
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10 R.L. Dewar & R.F. Abdullatif 

a slow drift in the phase while the dispersion dyadic d 2 u!k/d\s.dk causes 
spreading of the modulations. 

An amplitude-modulated wave can be represented as the sum of the 
unmodulated carrier wave and upper and lower sidebands 

A = Aq exp(— iAwrji) 

x {l + a + exp(iK-r - iilt) + a* exp(-iK-r + iQ*t)} , (15) 

where Aloq = Alu[\Aq\]. 

Linearizing in \a±\, \A\ = \Ao\[l + + a_) expi(K-r — Qt) + c.c], 
and using this in Eq. (|14(l we find 

= , (16) 
where 5wk (denoted aAw in Ref. is defined by 

Scjk = |-4olexp(— iK-r) f d 2 x^^- exp(iK> 
J o\A\ 







a+ 


\5lok 







r) . (17) 



Setting the determinant of the matrix in Eq. (|16(l to zero gives the 
dispersion relation for plane- wave modulations 

n - K -£)"-5 m &(** + 5 m &)- < 18 » 

The criterion for modulational instability is that O be complex, f2 = Q r +ir, 
r > 0, and from Eq. 118(1 we immediately see that this occurs, for sufficiently 
small K, if and only if there exist directions for K in which 

5. Nonlinear frequency shift 

The nonlinear frequency shift Aw in a general fluid or plasma is composed 
of two parts. The first is that due to the intrinsic nonlinearity of the medium 
and the second is that due to Doppler-like shifts [see Eqs. p0(l and (II ip ] 
associated with nonlinearly induced mean flows. 

However, in the case of drift or Rossby waves described by the CHM or 
MHM equations, the intrinsic nonlinear frequency shift is zero (or, at most, 
of higher order than quadratic). To see this, consider the terms in Eq. © de- 
scribing nonlinear wave- wave (including self) interactions: {tpi, tpi + S Sy ofpi} 
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and {"01, V 2, 0i}, where the Poisson bracket of two functions / and g is de- 
fined by 

{f,g} = zxV/'Vg = d x f.d y g - d x g.d y f . (20) 

Since we assume k y ^ 0, ip 1 is zero (to leading order, at least) and 
6:'i = As the Poisson bracket {/, /} = for any /, the first nonlin- 
ear self-interaction term vanishes. Similarly, because we are considering a 
monochromatic carrier wave, the second self-interaction term also vanishes 
to leading order: {tpi, V 2 ?/^} ~ — k 2 {ipi, ipx} = 0. 

For the calculation of the nonlinearly excited mean flows, we will need 
to evaluate the above term more accurately, which is best done via the 
useful identity 

{/, V 2 /} = d X dy [{d X ff - (dyff] - - C} 2 ) (dj.dj) . (21) 

(Some earlier discussion of this identity can be found in Ref. 1261 1 Averaging 
Eqs. Q and |JS} over the fluctuation scale, 

(d t + v*-V)(l - p 2 V 2 ) = p?P({&, V^i}) , (22) 
dt(ds,o - P 2 V 2 )(Jj) = P 2 P({^ V 2 ^}) . (23) 

For both the CHM and MHM cases, the small p 2 V 2 = 0(p 2 K 2 ) term 
on the LHS in Eq. (12211 is negligible compared with 1. 

5.1. Modulational instability for MHM equation 

However, on the LHS of Eq. l12.'i[l the leading term 1 does not occur in the 
MHM case, s = 1, so the p 2 V 2 term must be retained. Consequently, in this 
case (ip) is smaller than (ifi) by a factor 0(p 2 K 2 ). That is, the mean flow, 
oc zxK, is predominantly zonal so we lose no real generality in assuming 
K = Kx. and setting (ip) — 0. Then, using the identity Eq. (|21|l. dividing 
Eq. I|23|l by p 2 and integrating twice with respect to x we fino^ 

dtMx,t) = Pd^x.dy^t = 2k x k y \A\ 2 , (24) 

with the second form following from Eq. 

We can convert the time derivative to a spatial derivative by noting that 
the RHS of Eq. (jrl|) is small (assuming |A| is small) so, to leading order, 
dtA = — Vg-VA; that is, the modulations move at the group velocity. This 
also applies to quantities like i/jq driven by \A\, so, to leading order, Eq. (|24|) 
becomes 

dx ^ = -^y\A\ 2 . (25) 
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12 R.L. Dewar & R.F. Abdullatif 

By Eq. = d x iftoy, so Eqs. and (|12|) give the nonlinear frequency 

shift for case s = 1 as a simple function of \A\ 

Acu= {1+P l k2)2ky \A\ 2 . (26) 

As the nonlinear frequency shift in this case is a simple quadratic func- 
tion (rather than functional) of \A\, the modulated frequency shift parame- 
ter Slvk, Eq. I|17|) . in the nonlinear Schrodinger equation is now given sim- 
ply by Sujk = dAw/dln \A\ = 2Acl>o- The modulational instability criterion, 
Eq. ljT§|) . now becomes —k y d 2 ujk/dk 2 < 0. That is, from Eq. |JT3J|, the mod- 
ulational instability criterion for the modified Hasegawa-Mima equation 
case, s = 1, is 

1 - 3p 2 k 2 x + p 2 k 2 y >0, (27) 

which a gre es with Ref. and Ref. but not with Ref. OH who, due to a 
misprinP"! reproduced in Ref. EJ omit the factor 3 multiplying p 2 k 2 . (An 
apparently similar inequality to that in Ref. ED appears in Ref. but this 
is not really relevant as they consider only a drift wave propagating in the 
poloidal direction — their k x is our K x .) 

If criterion Eq. (|27|l is fulfilled, the growth rate curve, T 2 vs. K 2 , is an 
inverted parabola with maximum at 

P __ a K (l + P 2 k 2 f/ 2 \A\ 

max - AO; , A max - (i _ 3p2 ^ + ^2)1/2 p 2 v ^ ■ W 

This extends the small- K result in Eq. (15) of Ref. |S](who use the notation 
q for our K) to get a turnover in T at large K, as was also found using 
a mode-coupling approach by Chen et al® via the gyrokinetic equation in 
toroidal geometry within the ballooning approximation, and by Lashmore- 
Davies et using the modified Hasegawa-Mima equation. 

As the latter authors base their analysis on the same model as used in 
the present paper, we can make a precise comparison between our Eq. (|I8(I 
and their modulational dispersion relation, Eq. (43) in the small q = K x 
limit implied by our envelope approach. Expanding their quantities S± in 
q, it is easily seen that, to leading order, their 5+ + 6- = q 2 d 2 ujk/dk 2 
and (<5+ — S-)/2 = qduik/dk x , while their expression 2flQ\Ao\ 2 Fo(ko,q) is 
5u>k = 2Aojo, with Aluq given by Eq. 1)26(1 above. Completing the square in 
their Eq. (43), we see that correspondence between the two modulational 
dispersion relations is exact in the small q = K x limit. 
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Note also that the modulations and zonal flows have finite frequency, 
even without geodesic effectsf 621 as they propagate radially with phase ve- 
locity equal to the carrier group velocity, dujk/dk x = —Zp 2 k x k y v*/{\ + 
p 2 k 2 ) 2 from Eq. (Q2J. 



5.2. Modulational instability for CHM equation 

In the unmodified case, s = 0, there is actually no compelling reason to 
make the split into zonal and nonzonal components. Averaging Eq. 
over the fluctuation scale (noting that v^-VtA = and neglecting small 
terms) we find 

(d t + v*-V)Vo = Ps 2 (v E .VV 2 ^) = p 2 ({^, V 2 Vi}} ■ (29) 

Integrating from t = — oo, where ipo is assumed to vanish, along the trajec- 
tory of a fluid element moving at the drift speed and assuming the modula- 
tions in the forcing term on the RHS to be advecting at the group velocity 
we find 

M*, V, *) = Pl J drU^i, VVi}> (x - |gr, y + („, - |^)r, t 

X (30) 
Using the ansatz Eq. JSJ, the dispersion relation Eq. i|10[) . and the iden- 
tity Eq. (|21(l we find the nonlinear frequency shift to be the functional 



X 



Aw = T^T2 (*iA - W[(k 2 x - k 2 y )d x d y - k x k y (d 2 - d 2 )} 

Ps K 

1 



dr\A\ (x - v gx r, y + («* - v gy )r, t) . (31) 



Perturbing Aw with a small modulation 6\A\, replacing S\A\ by expiK-r, 
and substituting in Eq. (|17fl we get the modulated frequency shift 



5u) 



K 



^k 2 \A \ 2 (kyK x - k x K y )[(k 2 x - k 2 y )K x K y - k x k y (K 2 - K 2 ) 
l + p 2 k 2 K x v gx + K y (v gy - V*) 



Ug x T J-Iyyugy 

(32) 

Equation (|18|l then gives the dispersion relation for small modulations. 
Clearly, this is considerably more complicated than found in the MHM 
case and will not be analyzed further here except to make comparison with 
the results of Ref. 03 who take K y — 0. In this case 

5u, K = 2lM2K2x p 2 k 2 (l+p 2 k 2 ), (33) 
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which gives a modulational dispersion relation in the small K x limit in 
essential agreement with Eq. (19) of Ref.El who note that the modulational 
instability criterion is the same as that for the MHM, Eq. (|27ll . However, the 
resonance at K-(v g — v») = arising from the vanishing of the denominator 
in Eq. will give higher growth rates for oblique modulations, so it is 
not clear that this special case is of great significance for Rossby waves in 
the absence of boundaries. 

6. Conclusions 

We have derived a nonlinear Schrodinger equation for modulations on a 
train of drift or Rossby waves in a very universal, if heuristic, fashion. The 
nonlinear Schrodinger equation has been widely studied in other applica- 
tions and is known to have soliton solutions. However, we have analyzed 
it only for stability to small modulations and have found criteria in agree- 
ment with those found by Smolyakov et a/0 for modulation waves with 
zonal phase fronts. 

Our results are encouraging as a step towards explaining the experi- 
mental discovery by Shats and Solomon^ of modulational instability as- 
sociated with low- frequency zonal flows, but the Hasegawa-Mima equation 
is rather too simplified for direct comparison with experiment and further 
work remains to be done in this regard. 
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